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ABSTRACT 



We have mapped the central region of the Seyfert 1 galaxy NGC 1097 in 
12 CO(J = 2-1) with the Submillieter Array (SMA). The 12 CO(J = 2-1) map 
shows a central concentration and a surrounding ring, which coincide respectively 
with the Seyfert nucleus and a starburst ring. The line intensity peaks at the 
nucleus, whereas in a previously published 12 CO(J = 1-0) map the intensity 
peaks at the starburst ring. The molecular ring has an azimuthally averaged 
12 CO(J = 2 - 1)/(J = 1 - 0) intensity ratio {R21) of about unity, which is 
similar to those in nearby active star forming galaxies, suggesting that most of 
the molecular mass in the ring is involved in fueling the starburst. The molecular- 
gas-to-dynamical mass ratio in the starburst ring shows a somewhat lower value 
than that found in nearby star forming galaxies, suggesting that the high R21 of 
unity may be caused by additional effects, such as shocks induced by gas infall 
along the bar. The molecular gas can last for only about 1.2 x 10 8 years without 
further replenishment assuming a constant star formation rate and a perfect 
conversion of gas to stars. The velocity map shows that the central molecular 
gas is rotating with the molecular ring in the same direction, while its velocity 
gradient is much steeper than that of the ring. This velocity gradient of the 
central gas is similar to what is usually observed in some Seyfert 2 galaxies. To 
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view the active nucleus directly in the optical, the central molecular gas structure 
can either be a low-inclined disk or torus but not too low to be less massive than 
the mass of the host galaxy itself, be a highly-inclined thin disk or clumpy and 
thick torus, or be an inner part of the galactic disk. The R21 value of ~ 1.9 of 
the central molecular gas component, which is significantly higher than the value 
found at the molecular gas ring, indicates that the activity of the Seyfert nucleus 
may have a significant influence on the conditions of the molecular gas in the 
central component. 

Subject headings: galaxies: individual (NGC 1097), galaxies: active, galaxies: 
Seyfert, galaxies: nuclei, galaxies: ISM 



1. INTRODUCTION 



NGC 1097 (SB(s)b: lde Vaucouleurs et aDll99~lh is a nearby, inclin ed f46°;l0ndrechen. van der Hulst. fc 



19891 ). barred spiral galaxy at a distance of D = 14.5 Mpc (1"= 70 pc: lTullylll988l ). It hosts a 
Seyfert 1 nucleus, which is surrounded by a circumnuclear starburst r ing with a radius of 10" 
(0.7 kpc). NGC 1097 was originally identified as a LINER faeellll983h . but over the past two 
decades has shown Seyfert 1 activity as evide nced by the presence of broad double-peake d 



Balmer emission lines (FWHM « 7500 km s 1 : IStorchi-Bergmann. Baldwin, fc Wilson 



The starburst ring is rich in molecular gas (1.3 x 10 9 M Q ; iGerin et al.lll98 8). and exhibits a 
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hig h star formation rate of 5 yr~* as es timated from its extinction-corrected Ha luminos- 
ity (IHummel. van der Hulst. fc Keellll987l). A network of dusty spiral features connects th e 
starburst ring with the nucleus (IBarth et al.lll995l : iPrieto. Maciejewski. fc Reunanenl 120051 ). 
Optical spectroscopic imaging observations show evidence for radial streaming motion s as- 
sociated with the spiral structures leading to the unresolved nucleus (IFathi et al.l 120061 ). 



Kohno et al.l (120031 ) have previously mapped the central ~ V of NGC 1097 in 12 CO(J 
= 1-0) and HCN(J = 1-0) with the Nobeyama Millimeter Array (NMA). Their 12 CO(J 
= 1-0) map shows a ring-like structure consisting of several bright knots located at the 
starburst ring, and this molecular g as ring structure is similar to that revealed by the single 
dish 12 CO(J = 1-0) observations jGerin et al. l ll988h . The interferometric 12 CO(J = 1-0) 
map also revealed a relatively weak central concentration coincident with the nucleus. Their 
HCN(J = 1-0) map shows similar knots in the starburst ring, but in contrast to their 12 CO(J 
= 1-0) map a relatively bright central concentration of the HCN emission coincides with the 
nucleus. This implies that the molecular gas associated with the nucleus is relatively dense 
(riHj > 10 4 cm" 2 ). Simi lar den se central gas conc entrations have b een reported in N GC 



1068 flJackson et alJll993 l. M51 flKohno et al.lll996[ ). and NGC 6951 flKrips et al.ll2007f ). all 
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of which are type 2 Seyferts. In all these cases, the dense central gas concentrations have been 
attributed to a highly incline d circumnuclear molecular torus invoked by AGN unification 



models (e.g.. lAntonuccilll993l ) to obscure the central engine from direct view in the optical. 



This is the first time, however, that such a dense central gas concentration has been seen 
in a Seyfert 1 galaxy, warranting observations at higher angular resolutions to determine 
whether it is a good candidate for the hypothesized circumnuclear molecular torus. 

In this paper, we study th e central regio n of NGC 1097 in 12 CO(J = 2-1) using 
the Submillimeter Array (SMA; IHo et all 120041 ) . While the previous HCN(J = 1-0) and 
12 CO(J = 1-0) interferometric observations are suitable for tracing density variations in 
the molecular gas, they are not suitable for studying variations in the gas temperature. If 
the 12 CO(J = 1-0) line is optically thin, the intensity ratio between the 12 CO(J = 1-0) 
line and the higher- J l ines provides constraints on the gas temperature (and density) (e.g., 
Matsushita et al.ll2004l ). In addition, we observe the 12 CO(J = 2-1) line at a factor of 2 - 3 
higher angular resolution than the previous observations in 12 CO(J = 1-0) and HCN(J = 1 
- 0), allowing us to examine for the first time the spatial- kinematic structure of the central 
gas concentration. 



2. OBSERVATIONS AND DATA REDUCTION 

We observed the central region of NGC 1097 in the 12 CO(J = 2-1) line (rest frequency 
of 230.538 GHz) with the SMA, which has a primary beam of 52" (3.6 kpc) at this frequency. 
The receivers were turned to observe the 12 CO(J = 2-1) line at the upper side band. The 
observations were performed on 2004 July 23 and October 1 with eight 6-m antennas in 
the compact configuration. The 225 GHz zenith opacity of the two observations are ~0.15 
and ~0.3, and the T sySi DSB are ~ 200 K and ~ 350 K respectively. We placed the phase 
center at «2ooo = 02h46ml8 ! !96 and £2000 — -30°16'28'.'897, which corresponds to the position 



of th e AGN defined by the intensity peak of the 6-cm continuum emission (IHummel et al. 



19871 ). The SMA correlator has a total bandwidth of 2 GHz, and was configured to provide 
frequency resolution of 0.8125 MHz (~ 1 km s _1 ) and 3.25 MHz (4.2 km s" 1 ) for the July 23 
and October 1 observations respectively. We observed Uranus for bandpass and absolute flux 
calibrations. We used J0132-169 and J0423-013 for complex gain calibrations. J0132-169 is 
weaker but closer to the source (21? 5 away), and was used for phase calibration only; J0423- 
013 is stronger but further from the source (37° away), and used for amplitude calibration. 
The absolute flux accuracy and the positional accuracy have been estimated as ~15% and 
O'.'l, respectively, in these observations. 



The observation on July 23 was binned to the same velocity resolutions as that on 
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October 1 (3.25 MHz), and then combined with each other. The data were calibrated 
with the Owens Valley Radio Observatory software package MIR, which was modified for 
the SMA. The images were CLEANed using the NRAO software package AIPS, and with 
natural weighting have an angular resolution of 4"1 x 3"1 (290 pc x 220 pc) at a position 
angle (P. A.) of 168°. The rms noise level of the individual channel maps is 75 mJy beam"" 1 . 
We used the "MOMNT" task in AIPS to make the integrated intensity and the intensity- 
weighted mean-velocity maps. The line free channels were binned to create a continuum 
map with a rms noise level of 4.4 mJy beam -1 ; no significant emission was detected in this 
map. 



RESULTS 



We detected the 12 CO(J = 2-1) line spanning a total line width of ~ 570 km s -1 
at the 3a detection threshold in the ch annel maps. The line w idth is comparable to that 
measured in the single-dish observation (IPetitpas fc Wilsonll2003l ) of ~ 550 km s -1 using the 
JCMT, which has a primary beam of 21", although the integrated line intensity in our map 
is only 60 ± 13% of that measured in this single-dish observation (see Sect. 13.31) including 
uncertainties of the flux in both observations. 



3.1. Spatial Distribution 



We show the 12 CO(J = 2-1) integrated intensity map in Figure d](a). A relatively 
strong molecular concentration (hereafter, central component) is detected at the very center 
in a region of a diameter of 350 pc (5"), surrounded by a weaker molecular ring of a diameter 
of 1.4 kpc (20"). The centroid of the central molecular gas component is shifted by ~ 0"7 
from t hat of the radio cont inuum core at 6 cm, which presumably marks the location of the 
AGN (IHummel et al.lll987t ). We shall return to this positional difference in Section [3T2l 



The molecular ring coincides with the starburst r i ng, c onsistent with the previous 
12 CO(J =1-0) maps JGerin et al l Il988l: iKohno et alj liooih. Unli ke the tightly-wound 



spiral structure seen in optical maps (e.g., iRickard 



1975 



Barth et al. 



1995L the molecular 



ring traces a complete circle as seen in the J — K$ color map of iPrieto et al.l (120051 ) . The 
molecu lar ring is composed o f several knots that closely resembles the 1.5 GHz continuum 



map of IHummel et al.l (119871 ). In addition, we detected weak molecular emission extending 
from the NE and SW of the molecular ring, coinciding with dust lanes along the large stellar 
bar (indicated by the two straight lines in Fig. (T^,). 
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In Figure [T^b), we show the 12 CO(J = 1-0) map made by lKohno et al.l (120031 ) with the 
NMA. As can be seen, the molecular ring is more clearly discernible and also more clearly 
separated in our 12 CO(J = 2-1) map than in the 12 CO(J = 1-0) map. The intensity peak 
in the 12 CO(J = 2-1) map is located at the nucleus, whereas that in the 12 CO(J = 1-0) 
map is located at the NE part of the ring. The SW side of the ring is brighter than the NE 
side in the 12 CO(J = 2-1) map, but the reverse is true in the 12 CO(J = 1-0) map. Both 
of the molecular gas peaks at the NE and SW sides of the ring are located at the regions 
where the dust lanes connect with the ring, and these molecular gas peaks are often seen 
in ba rred spiral galaxies, which are known as the twin-peak morphology (e.g.. iKenney et al. 



19921). 



3.2. Kinematics 



In Figure [2(a), we show the intensity-weighted 12 CO(J = 2-1) mean-velocity map. 
Both the central component and molecular ring show an overall velocity gradient in the NE 
to SE direction alo ng a P.A. of 135°, that is similar to the major kinematic axis of the large- 
scale galactic disk (jOndrechen et al.lll989l ). The emission is blueshifted on the northwestern 
side of center, and redshifted on the so utheastern side of c enter with respect to the systemic 
velocity of 1254 km s _1 for NGC 1097 (jKohno et al.ll2003l ). The gas motion in the molecular 
ring appears to be dominated by circular motion (i.e., isovelocity contours pe r pendi cular 
to major axis), whereas in the 12 CO(J = 1-0) velocity map of IKohno et al.l (120031 ) the 
isovelocity contours have a symmetric S-shape with the end of the S-shape nearly parallel to 
the dust lanes along the large stellar bar. This non-circular motion is not as prominent in 
our 12 CO(J = 2-1) map perhaps because the aforementioned dust lanes are not as strongly 
detected in 12 CO(J = 2-1), along with the fact that th ey are closer to the edge of our 
primary beam than in the 12 CO(J = 1-0) observations of IKohno et al.l (12003 ). 

A position- velocity diagram (hereafter p-v diagram) of the 12 CO(J = 2-1) emission 
along its major kinematical axis (P.A. = 135°) is shown in Figure[2](b). Positive and negative 
velocities correspond to redshifted and blueshifted velocities, respectively, relative to the 
systemic velocity. The rapidly rising part of the rotation curve corresponds to the central 
component, and the flat part to the molecular ring. The rotation curve is symmetric on both 
sides of the center, rising steeply to ±235 km s _1 at ±2'/5 (±175 pc) and flattening outside 
±2'/5. This indicates that the size of the central component is about 350 pc in diameter. 
As can be seen in Figure 12(b), the emission from the central component is stronger on the 
redshifted southeastern part, causing the centroid of the emission to be shifted by CP.7 towards 
the southeast as mention in Section I3TT1 (and shown in Fig. [Ih). The velocity gradient of the 
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central component is 1.31 ± 0.14 km s 1 pc 1 (Fig. [Tb), which is similar with that derived 
in th e optical from ionized gas (~ 1.1 km s _1 pc -1 ; IStorchi-Bergmann. Wilson. &: Baldwin 
1996h . 



3.3. 12 CO(J = 2 - 1)/(J = 1 - 0) Line Ratios 



Our 12 CO(J = 2 - 1 ) map detected about 60% of the flux measured with the JCMT 
dPetitpas fc Wilson|[2003h . The shortest projected baseline in our SMA observations is 4.8 kA 
(~ 6 m), and so the largest structure we can detect is ~52", comparable in size to the SMA 
primary beam. We therefore assume that the molecular gas resolved out in our map is 
uniformly distributed in space and in velocity over the entire line width, and have corrected 
for ~ 40% of the missing flux. We applied a primary beam correction to both the 12 CO(J 
= 2-1) and 12 CO(J = 1-0) maps, and convolved the derived moment maps to the same 
angular resolution of 6"5 x 3"5 with P. A. of 0°. These moment maps were then used to derive 
the 12 CO(J = 2 - l)/ 12 CO(J = 1-0) line ratios. Note that we also calculated the line ratios 
with itu-matched data, and the ratios are consitent with the missing flux corrected line ratios 
within errors. To derive the molecular gas column densities and masses, we need to use the 
missing flux corrected data (see Sect. I3.4p . To make consistency within this paper, we use 
the missing flux corrected data for the following calculations. 

In Figure [3], we show the azimuthally averaged radial intensity distributions of the 
12 CO(J = 2-1) (dashed line) and 12 CO(J = 1-0) (dotted line) lines, and the azimuthally 
averaged 12 CO(J = 2 - l)/ 12 CO(J = 1-0) line intensity ratio (.Rai) (solid line). The radial 
intensity profile in both 12 CO(J = 2-1) and 12 CO(J = 1-0) peak at the center, correspond- 
ing to the central component, and exhibit a secondary peak at a radius of 10", corresponding 
to the molecular ring. The line ratio R21 exhibits a similar behavior, peaking at the central 
component and exhibiting a secondary peak at the molecular ring. At the nucleus, R21 is 
derived as 1.9±0.2 at a beam size of 6"5 x 3'.' 5, which corresponds to the approximate angular 
extent of the central component along its major kinematic axis. In the molecular ring, the 
azimuthally averaged value between radii of 8" and 12" is R21 = 1.3 ± 0.2. 

The averaged R21 at the molec ular ring is similar to or somewhat larger than the global 
i n spiral galaxies of ~ 0.7 — 0.9 (IBraine et al.lll993l ; lLavezzi et al.lll999l ; lHafok fc Stutzki 
20031). The -R21 at gia nt molecular clouds (GMCs) in the nearby star forming r egion Orion 



(Sak amoto et al.lll994j ) or GMC-scale R21 in the nearby starburst galaxy M82 (IWeiss et al. 
200 ll ) also shows similar R21 of about unity. On the other hand, the R21 at the central com- 



ponent is about twice as high as those in the abovem entioned sources, and consistent with the 



ratio reported in the Seyfert 2 galaxy NGC 1068 (IBaker fc Scovilld Il998l ; ISchinnerer et al. 
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2000l ). The R21 at the central component is also similar to that found at active star form- 



ing regions or at the interfaces of molecular clouds and ionized gas in the near by barred 
spiral galaxy IC 342 ([Turner. Hurt. &: Hudson! Il993l ; iMeier. Turner, fc Hurtll20001). Galactic 



objects, such as mole cular outflows from young stellar objects (e.g. 



Richardson et al 



19851: 



Chandler et afll996h or molecular gas around supernova remnants (Ivan Dishoeck. Jansen. fc Phillips! 



1993 



Seta et al.lll998l ) also show R21 > 2. Note however that the high value of R21 are not al- 



w ays seen in active galaxies; for instanc e, interferometric observations of NGC 3227 (Seyfer t 



1; I Schinnerer. Eckart. fc 



and NGC 6574 (Seyfert 2; 



accom 



200oh. NGC 37 18 (LINER/Seyfert 1: iKrips et al.l bpQSh . 



Lindt-Krieg et al.ll2008l ) show R21 of around unity. 



3.4. Physical Properties of the Molecular Gas 



We derived the column density of the molecular gas in the central component and 
molecular ring assuming the standard Galactic conversion factor of 3.0 x 10 20 cm -2 (K km 
s -1 ) " 1 between the 12 CO(J = 1-0) line inten sity and column density of molecular hydrogen 
gas (jScoville et al.l 1 198 71 ; ISolomon et al.lll987l ). The 12 CO(J = 2-1) integrated intensity at 
the center (beam size of 3"1 x 4"1) is 376.7±7.2 K km s _1 . Assuming the same line ratios 
derived from a beam size of 6'.'5 x 3"5, namely, R21 = 1.9 ± 0.2 for the central component, 
we derive a column density for this component of (5.9 ± 0.6) x 10 22 cm~ 2 , and the molecular 
hydrogen gas mass of (6.5 ± 0.7) x 10 7 M . In the molecular ring, the azimuthally averaged 
12 CO(J = 2-1) integrated intensity between the radii of 8" and 12" is 135.7 ± 7.2 K km 
s -1 . With R21 = 1.3 ±0.2, the azimuthally averaged column density at the molecular ring is 
therefore (3.0±0.3) x 10 22 cm -2 , and the molecular hydrogen gas mass is (5.8±0.6) x 10 8 M . 
Here we compare our estimated molecular gas mass of the ring with that estimat ed using the 



single- dish data. We re-calculated the gas mass inside 20" in radius derived by iGerin et al. 



( 119881 ) using the conversion factor mentioned above, and it is calculated as 1.1 x 10 9 M . We 
then subtract the gas mass of the center from this gas mass. The gas mass of the center is 
ca lculated as 2 . 4 x 10 8 M using the central position data with a of ~32 K km s -1 (Fig. 3 



of Gerin et al. ( 19881 )). The gas mass of the ring is therefore calculated as ~8.6 x 10^ 



M 



©• 

This value is somewhat larger than ours, and this may be due to their larger radius of 20", 
and therefore their results may be detecting emission outside the ring. 

To determine the tem perature and density of the molecular gas, we use the LVG method 
( iGoldreich fc Kwanl Il974l ). We assume a one zone model, which assumes that both the 
12 CO(J = 2-1) and 12 CO(J = 1 - 0) emissions orig i nate from the same region. The 
collision rates for CO are taken from iFlower fc Launayl (119851 ) for the temperature regime 
T = 10 - 250 K, and iMckee et all (|l982j) for T = 500 - 1000 K. We first assume a stan- 



- 8- 



dard Z( 12 CO)/(dv/dr) of 5.0 x 1(T 5 (km s" 1 pc" 1 )- 1 , where Z( 12 CO) = [ 12 CO]/[H 2 ] is the 
abundance ratio and dv/dr is the velocity gradient of the molecular gas. 

We then compute from the LVG model the line ratio R21 and 12 CO(J = 1-0) opacity as 
a function of molecular hydrogen number density nn 2 and kinetic temperature Tk as shown 
in Figure|H(a). The result for R21 = 1.9 ±0.2 as inferred for the central component indicates 
that T k > 400 K, and n H2 ~ 3 x 10 4±1 cm" 3 , and 12 CO(J = 1-0) opacity below unity 
The density estima tion is consistent w ith the detection of HCN(J =1-0) line from the 



central component (jKohno et al.ll2003l ). which indicates the molecular gas density as high as 
riH 2 ~ 10 4 cm -3 . Note that the estimated kinetic temperature is highly dependent on the 
assumed Z( 12 CO)/ (dv/dr), namely highly dependent on the [ 12 CO]/[H 2 ] relative abundance, 
on the velocity gradient, or on both. Fixing tzh 2 = 1 x 10 4 cm -3 , we plot R21 as a function of 
Tk and Z( 12 CO) /(dv/dr) in Figure IU(b). As can be seen, R21 increases roughly linearly with 
Z( 12 CO) /(dv/dr), and around the standard Z( 12 CO) /(dv/dr) of ^10~ 5 (km s" 1 pc" 1 )" 1 , 
we find that a kinetic temperature of at least 100 K is required to reach f? 21 of 1.9 ± 0.2. If 
on the other hand Z( 12 CO) /(dv/dr) is an order of magnitude lower than the standard value 
(i.e., order of 10~ 6 ), the kinetic temperature would be in the range of Tk ~ 30 — 250 K, which 
is comparable with the temperature range normally found in molecular clouds. 

For the molecular ring, which has an average R21 = 1.3 ± 0.2, Tk > 100 K and tih 2 ~ 
1 x 10 4±1 cm -3 assuming the standard Z( 12 CO) /(dv/dr). Again, these values are sensitive 
to the assumed Z( 12 CO)/ (dv/dr), and if this value decreases, T k and n Hz also decrease. 



4. DISCUSSIONS 

4.1. Molecular Gas in the Starburst Ring 

The average i? 2 i of about unity in the molecular ring suggests that the overall properties 
of the molecular gas in the starburst ring is optically thick and thermalized. As mentioned 
in Section 13.31 this value is similar to that of the global R21 of other spiral galaxies or that in 
the active star forming GMCs. Similarity of the molecular gas ratios between the molecular 
gas in the ring and the molecular clouds at active star forming regions indicates that the 
molecular gas in the ring is closely related to the starburst activities occurrin g in the ring. 



This is also supported by the high HCN(J = 1 - 0)/ 12 CO(J = 1 - 0) of 0.16 dKohno et al. 



20031 ). Since stars form from molecular gas, we believe that this molecular gas is actually 
fueling the starburst activities in the ring. The gas mass content in the starburst ring is, 
on the other hand, not as high as star forming galaxies, about 8% of the dynamical mass; 
the dynamical mass within the radius of 8" - 12" (560 - 840 pc; roughly corresponds to the 
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width of the molecular ring) can be estimated as (6.9 ± 0.2) x 10 9 M Q , using the rotational 
velocity of 235 km s -1 (Sect. 13.21) and as suming the inclination of the molecular ring is the 
same as that of the host galaxy of 46° (jOndrechen et al.lll989l ). The molecular gas mass 
in the ring is 5.8 x 10 8 M (Sect. I3.4p . so that the molecular gas mass to dynamical mass 
ratio is calculated as (8 ± 1)%. This value is a bit l ower than the th r eshol d to have star 
formation activities at nuclei of ~ 10% suggested by ISakamoto et al.l (119991 ). One reason 
for the discrepancy between high R 2 \ and HCN(J = 1 - 0)/ 12 CO(J = 1 - 0) rat ios and low 



molecular gas mass to dynamical mass ratio may be that ISakamoto et al.l (119991 ) are taking 
mass ratios within 500 pc in radius, but we are taking the mass ratio in a certain range of 
radius. Another possibilities is that the gas in the staburst ring is highly affected by shock 
induced by the gas infall along the bar, so that the high intensity ratios may not tracing the 
star formation activities, but tracing shocked gas. 

The upper li mit of the mass of n ewly formed (1.5 Gyr old) stars in the star-forming 
ring is ~ 10 9 M Q (jQuillen et al.lll995l ). which is roughly comparable to or somewhat larger 
than the mass of molecular gas (5.8 x 10 8 M Q ). This implies that a large fraction of the 
molecular gas has already turned into stars within a timescale of 1.5 Gyr. We can estimate 
the molecular gas consumption timescale with considering the star formation rate (SFR) 
and t he molecular gas ma ss of the ring. The SFR in the molecular ring is about 5 M 
yr _1 (jHummel et al.lll987l ). Without further replenishment, the molecular gas in the ring 
can last for only about (1.2 ± 0.1) x 10 s years assuming a constant SFR and a perfect 
conversion of gas to sta r s. Th is value is roughly an order of magnitude lower than the value 
derived by iGerin et al.l (119881 ). which is due to their somewhat larger molecular gas mass, 
inclusion of atomic hydrogen gas mass assuming it has the same amount as the molecular 
gas mass, and using lower SFR of ~ 2 M yr _1 . In addition to the gas consumption by 
star formation, the amount of molecular gas in the ring also decreases by the infall toward 
the nucleus. As mentioned in Section [Tj, the gas infalling motion has been observed from 
the ring to the nucleus. The value for the mass infall rate is not known, but this effect 
obviously decreases the lifetime of the molecular ring. On the other hand, gas replenishment 
by the infall along the large scale bar from the outer part of the galaxy to the molecular ring 
increases the lifetime. The value for the mass infall rate is again not known, so the lifetime 
of the molecular ring is decided by the balance between these effects. 



4.2. Nature of the Nuclear Component 



Our 12 CO(J = 2-1) results show that the molecular gas located at the center is spatially 
coincident with the AGN, which is consistent with the previous interferometric maps of the 
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12 CO (J = 1 - 0) and HCN(J = 1-0) lines flKohno et al.l 120031 ). In addition, the central 
component is kinematically symmetric at the center with steep velocity gradient, suggesting 
that the molecular gas is rotating around the AGN. Here we discuss the relation between the 
central component and the AGN activities based on the line ratios of the central component 
and its physical conditions, which were derived using the LVG analysis results presented in 
the previous section. 

The 12 CO(J = 2-1) map peaks at the nucleus, different from the 12 CO(J = 1-0) 
map, and the R21 shows 1.9 ± 0.2. This value is significantly higher value than the global 
R21 of spiral galaxies or R21 in star forming GMCs in Orion or M82, but similar to the R21 
of molecular gas at jets, at supernova remnants, or at the surface between molecular gas 
and ionized gas (Sect. 13.31) . These results suggest that the high R 2 \ in molecular gas can be 
related to irradiation of UV photons from star forming regions, shock caused by interaction 
between molecular gas and outflowing or expanding materials, but not to the global galactic 
characteristics or activities. Hence similar activities can be the cause of high R 2 ± in the 
central region of NGC 1097. However, other activities that cannot be seen in star forming 
galaxies or in our Galaxy can also be the cause of high R^ , such as strong X-ray radiation 



from the Seyfert 1 nucleus fe.g..lLepp fc Dalgarnolll996l ; iKohno et al.ll200ll ; lUsero et al.ll2004 



Kohnoll2005l : I Meiierink. Spaans. fc Israeli l2007h . 



From LVG analysis, high temperature (> 400 K) and high density (~3x 10 4±1 cm~ 3 ) 
conditions are derived for the central component, assuming a standard Z( 12 CO)/(dv/dr). 
As ment ioned in Section I3"T4"| the density is consistent with the centrally peaked HCN(1 - 
0) map ( Kohno et al. 20031 ). Furthermore, these conditions are supported by results from 
infrared observations: Strong molecular hydrogen line, H 2 1 - S(l), is detected toward the 
nucleus without any detection of Br7 line, and the JHK images show red colors toward the 
nucleus, suggesting tha t the presence of hot dust coexisting with the dense molecular gas 
(IKotilainen et al.l 120001 ) . Such high density and high temperature conditions for molecular 
gas arou nd Seyfert nuclei are al s o der ived in other gal axies, such as the Se y fert 2 nucleus 



of M 51 (Matsushita et all Il998l . |200J) or NGC 1068 flRotaciuc et al.l 1 199 it iTacconi et al. 
99 4]), but not for the molecular gas in non-active galaxies such as IC 342 or our Galaxy 



([Matsushita et al.lll998l ). These physical values are, however, sensitive to Z ( CO) / (dv / dr 



namely to the 12 CO abundances, to the velocity gradient, or to both. As shown in Section [ 
lower Z( l2 CO)/(dv/dr) makes the derived temperature lower. This is because the opacity 
is linearly related to Z( 12 CO)/(di>/dr), and it will decrease if Z( 12 CO)/(dv/dr) decreases. 
Under low Z( 12 CO)/(dt>/dr) conditions, if the temperature rises, the 12 CO(J = 1-0) line 
can easily be optically thin, and easier to have high R 2 \ at lower temperature than the 
normal Z( 12 CO)/(dv/dr) condition. In either case, the molecular gas around Seyfert nuclei 
seems to have different properties from other normal environment, which seems to be largely 
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related to the Seyfert activities. Combined with the spatial and kinematical information, 
we suggest that the central component is kinematically and physically related to the Seyfert 
nucleus. 



4.3. Is the Circumnuclear Gas the Hypothetical Circumnuclear Torus? 



Given that the central component is closely related to the AGN, is this then the circum- 
nuclear molecular torus predicted by AGN unified models? Since NGC 1097 hosts a Seyfert 
1 AGN, we expect that the rotating circumnuclear gas is in a nearly face-on configuration 
to provide an essentially unobscured view to the bro ad-line region (BLR), if we apply the 
general unified model for AGNs (e.g.. lAntonuccilll993l ). Kinematics of the central component 
we observed are, h owever, similar to what is usually observed in other Seyfert 2 galaxies , 



such as NGC 1068 (I Jackson et al.lll993l : ISchinnerer et all 120001 ) or M51 (IKohno et al.lll996 



Scoville et al.lll998l ); the central component shows steep velocity gradient that can be ex- 
plained by edge-on disk or torus rotating around the Seyfert nucleus. In addition, we derived 
a molecular hydrogen column density (A% 2 ) toward the nucleus of (5.9 ± 0.6) x 10 22 cm -2 
(Sect. I3.4p . or 1.2 x 10 23 cm -2 in atomic hydrogen column density (Ah), which is about 
two orders of magni t udes larger than Ah derived from X-ray spectra of 1.3 x 10 
(jlyomoto et al.lll996l ; iTerashima et al.l 120021 ). 



Jl cm 2 



From these "inconsistent" results between our molecular gas observations and other 
observations at different wavelengths/frequencies, the structure of the central component 
can have two basic configurations; one is a nearly face-on disk-, ring-, or torus-like structure 
with a very fast rotation velocity, and another is a nearly edge-on disk-, ring-, or torus- 
like structure with a thin or clumpy structure. Due to our large synthesized beam size, 
we could not distinguish these two possibilities observationally, and also could not evaluate 
the thickness of the structure. Here we discuss the advantages and disadvantages of both 
possibilities. 

The former molecular gas configuration has a nearly face-on structure, so the direct view 
to the BLR is secured. But since the observed rotational velocity width has 470 km s -1 , 
the inclination corrected rotational velocity width for this configuration has to be 470 sin(i), 
where % is the inclination (0° corresponds to a face-on configuration). The nearly face-on 
configuration therefore leads to a high rotational velocity (e.g., 940 km s _1 even for the 
inclination angle of 30°). Note that this rotation is rigid r otation, and it is rare to see rigid 



Rubin et al. 


1982; 


Sofue et al. 


1999) 



in other galaxies (e.g., iRubin. Thonnard. &: Fordlll980 



and 5°, the dynamical mass within the central gas, namely within the radius of 2'.'5 (175 pc) 
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can be estimated to be 8.8 x 10 9 M , 7.3 x 10 10 M , and 2.9 x 10 11 M , respectively. The 
total dynamical mass of NGC 1097 within a radius of 7!5 (31.5 kpc) estimated from th e 



large-scale atomic hydrogen observations is (5.0 ± 0.8) x 10 11 M Q (IHigdon &: Wallinl 120031 ) . 
A disk with the inclination angle of 5° is impossible, since the dynamical mass of the central 
gas is almost the same with the mass of the whole galaxies. A disk with the inclination angle 
of 10° is still too large, since more than one tenth of the total mass is concentrated within a 
thousandth of radius. A disk with the inclination angle of 30° can be possible. We therefore 
think that this face-on configuration can be possible, only if the inclination angle is ~ 30° 
or larger. 

The latter molecular gas configuration has an edge-on structure, so the structure should 
have thin disk- or ring-like structure, or it can be torus-like structure but has to have a clumpy 
internal struc ture to expose the BL R at the center. Such clumpy structure is suggested 



theoretically (IWada fc Norman! 120021 ) . and with their model, the column density of ~ 10 21 
cm -2 is possible even with the inclination angle of ~ 60°. Under this model, the difference of 
the column density derived from our CO data and from the X-ray data can be explained; the 
spatial resolution of our observation is about 250 pc, which smears all the internal clumpy 
structures, and therefore the derived column density will be the average value and higher 
than that derived from X-ray observations, which only trace very narrow column densities 
due to the very small size scale of the X-ray emitting region (< 1 pc). Under these nearly 
edge-on configurations, the rotational velocity is in the typical values for other galaxies 
(< 500 km s _1 ), and therefore we do not need to invoke any special conditions. 

Of course, there is another possibility that the central component is nothing related 
to the hypothetical torus, namely the central component we observed is just a part of the 
galactic disk gas. This is supported by the similar trend and the smooth connection of the 
molecular gas kinematics at the nucleus and the ring. Since the infalling motion exists from 



the ring to the nucleus along the nuclear bar or spiral (jPrieto et al.ll2005l ; iFathi et al.ll2006l ). 
it is natural to pile up the gas around the nucleus with similar kinematics as that of the 
ring. In this case, gas will rotate around the nucleus with similar inclination as the outer 
disk or the ring, and it is natural not to cover the line-of-sight toward the AGN, namely the 
gas configuration will be similar to the Seyfert 1 nucleus. In addition, if the molecular gas 
piles up at the diameter smaller than our beam size of 200-300 pc, the discrepancy of the 
column density estimated from our data and the X-ray data can also be explained. Here we 
briefly estimate the necessary gas infall rate to create the central component with gas infall 
from the molecular ring. The H 2 masses of the central component is 6.5 x 10 7 M , so that 
the mass infall rate of 0.5 M yr _1 is required to create this component within the molecular 
ring gas consumption timescale of 1.2 x 10 8 years (Sect. 14.11) . Note that the line ratio of the 
central component is obviously different from that in other regions or other galaxies, and 
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this can be explained either by the shock caused by the gas inflow along the nuclear bar 
or nuclear spiral, or by AGN activities, such as irradiation of strong X-ray emission to the 
central component. 

In summary, the central component can be the hypothetical circumnuclear disk or torus 
with nearly edge-on (> 60°) clumpy structure, or less likely nearly face-on (> 30°) disk/torus. 
The central component can also be nothing related to the hypothetical disk/torus, and 
possibly created by the gas inflow from the molecular ring. The physical conditions of the 
central component are different from that of molecular gas in other regions, so that even it 
is not the hypothetical disk/torus, it should be highly related to the AGN activities or gas 
inflow toward the nucleus. 



5. SUMMARY 

We successfully imaged the central component and the molecular ring, which is spatially 
coincident with the AGN and the starburst ring, toward the central region of the Seyfert 1 
galaxy NGC 1097 using the SMA in the 12 CO(J = 2-1) line. Here are the summary for 
the nature of the central component we observed: 

1. We found that the 12 CO(J = 2-1) map shows an intensity peak at the central com- 
ponent, and different from the 12 CO(J = 1-0) map, which shows the intensity peak 
at the molecular ring. 

2. The 12 CO(J = 2 - l)/ 12 CO(J = 1-0) line intensity ratio for the central component 
is 1.9 ± 0.2, which is different from the global values in GMCs or in spiral galaxies of 
about unity or less. From the LVG analysis, we estimated that the central component 
is warmer (T K > 400 K) and denser (nu 2 ~ 3 x 10 4 cm -3 ) than that of the normal 
molecular clouds assuming a normal Z( 12 CO)/ (dv /dr). These values depend highly on 
the Z( 12 CO)/(di>/dr), and lower Z( 12 CO)/(di>/dr) for about an order of magnitude 
decreases the temperature and density of about an order of magnitude. The effect of 
intense star formation and/or AGN activities, such as shocks induced by numerous 
supernova explosions or strong X-ray radiation from AGN are presumably the causes 
of the unusual line ratio. 

3. Faster rotation feature of the central gas is observed in NGC 1097, which is similar 
results as observed in other Seyfert 2 galaxies. We interpret this feature as a highly- 
inclined clumpy disk/torus or thin disk to explain the fast rotation and the difference 
between the column density derived from CO and X-ray observations. A low-inclined 
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(i ~ 30°) thick disk is possible, but lower inclination than this value is less likely since 
it is rare to see a velocity > 500 kms -1 in nearby galaxies, and the total mass of central 
disk turns to be too large compared with the total mass of the galaxy. On the other 
hand, the central component can also be interpret as just an inner extention of the 
galactic disk, possibly created by the gas inflow from the molecular ring. 

4. Combining kinematical and spatial information with the physical conditions, we sug- 
gest that the central gas is related to the Seyfert activities or gas inflow. 

Here are the summary for the nature of the molecular ring we observed: 

1. The R21 of the molecular ring is 1.3 ± 0.2, which shows a similar properties to the 
star forming GMCs in our Galaxy or nearby starburst galaxies. In addition, since the 
molecular ring shows a good spatial coincident with the starburst ring, so we expect 
that the molecular ring is actually fueling the starburst activities. The molecular 
gas mass content with respective to the dynamical mass, on the other hand, shows 
somewhat low value, so that high R 2 i may not be related to star formation activities, 
but to other activities, such as shock induced by gas infall along the bar. Without 
further replenishment, it can last for only about 1.2 x 10 8 years, but since there seems 
to have gas flows from the dust lane along the large-scale bar toward the molecular ring, 
and from the molecular ring toward the nucleus, this timescale is highly uncertain. 
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comments. We also thank for the SMA staff for maintaining the operation of the array. The 
Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory 
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Fig. 1. — (a) 12 CO(J = 2-1) integrated intensity (moment 0) map of the central region of 
NGC 1097. The contour levels are (2.5, 5, 10, 20, 30, 40, and 50) x 3.7 Jy beam" 1 km s _1 , 
where la is 3.7 Jy beam -1 km s -1 . The synthesized beam is 3"1 x 4'.'1 with a P.A. of 168°, 
which is shown in the bottom-left corner of the map. The large white circle is the primary 
beam size (52") of the SMA. (b) 12 CO(J = 1-0) integrated intensity map. The contour 
levels are (2.5, 5, 7, 9, and 11) x 3.9 Jy beam -1 km s -1 , where la is 3.9 Jy beam -1 km 
s . The synthesized beam is 6"0 x 3"0 with a P.A. of 0°, which is shown in the bottom-left 
corner. The large white circle is the primary beam size (63 ") of the NMA . In b oth images, 
the central crosses mark the peak of the 6 cm continuum (IHummel et all 119871 ). the small 
red circles at a radius of 10" mark the radius of the molecular ring, and the two red lines 
mark the positions of the dust lanes. 
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Fig. 2. — (a) Intensity-weighted mean velocity (moment 1) map of 12 CO(J = 2-1). Large 
and small circles are marked as the primary beam size of the SMA (52") an d the size of 



starb urst ring (radius of 10"), respectively. The systemic velocity (1254 km s . iKohno et al. 



20031 ) is marked. The contour interval is 15 km s 1 . Dotted and solid curves are the 



blueshifted and redshifted velocities, respectively. The solid lin e labeled as "A" rep resents 
the direction of the major axis (P. A. is 135°) of NGC 1097 jHummel et alJll9R7h . The 
nuclear gas is rotating with the molecular ring in the same sense, while it is rotating faster 
than the ring, (b) Posit ion- velocity diagram of the 12 CO(J = 2-1) emission cut along the 
major axis (represented as "A" in the above figure). Contour levels are 3, 5, 10, 15, and 20a, 
where la = 75 mJy beam -1 . The zero velocity corresponds to the systemic velocity of 1254 
km s _1 marked as horizontal solid lines. The offset of 0" corresponds to the galactic center, 
and the LSR velocity of km s _1 corresponds to the systemic velocity. The fitted velocity 
gradient (1.31 km s -1 ) of the nuclear component is marked as a dash line. Note that the 
inclination of the galactic disk is not corrected. 
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Fig. 3. — Azimuthally averaged radial intensity distributions of the 12 CO(J = 2-1) (dashed 
line) and 12 CO(J = 1-0) (dotted line) lines in brightness temperature scale. The temper- 
ature scale is shown in the left-hand side of the vertical axis. Azimuthally averaged radial 
distribution of i? 2 i is marked as a solid line. The ratio scale is shown in the right-hand side 
of the vertical axis. 
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Fig. 4. — (a) i?2i is plotted as a function of kinetic temperature and H 2 number density 
riH 2 & t a standard Z( 12 CO)/(dv/dr) of 5.0 x 10~ 5 (km s _1 pc -1 ) -1 . Solid lines are the R21 of 
the nuclear component (1.9±0.2) and the molecular ring (1.3±0.2). Opacity for 12 CO(J = 
1 - 0) is marked as dashed lines, (b) R21 of the nuclear component is plotted as a function 
of T k and Z( 12 CO)/(dv/dr). Solid lines are the observed R 2 \ of the nuclear component with 
nn 2 fixed at a value of 1.0 x 10 4 cm~ 3 . 



